Creep of superconducting vortices in the limit of vanishing temperature: A 
fingerprint of off-equilibrium dynamics. 



o : 
o . 

(N ■ 
>v 



m 
> 
oo 
in 

(N 

O 
O 
O 



c 

o 
o 



13 



Mario Nicodemi a,h and Henrik Jeldtoft Jensen a 
a Department of Mathematics, Imperial College, 180 Queen's Gate, London SW7 2BZ, UK 
b Universitd di Napoli "Federico II", Dip. Scienze Fisiche, INFM and INFN, Via Cintia, 80126 Napoli, Italy 

We theoretically study the creep of vortex matter in superconductors. The low temperatures 
experimental phenomenology, previously interpreted in terms of "quantum tunnelling" of vortices, 
is reproduced by Monte Carlo simulations of a purely "classical" vortex model. We demonstrate that 
a non-zero creep rate in the limit of vanishing temperature is to be expected in systems with slow 
relaxations as a consequence of their off-equilibrium evolution in a complex free energy landscape. 
PACS numbers: 74.60.Ge 47.32.Cc 74.50.+r 75.45.+j 



There exist an abundance of experimental evidence 
that the relaxation rate of the magnetisation in type II 
superconductors (ranging from conventional, to organic, 
heavy fermions and a variety of H — T c superconductors 
does not vanish as the temperature T is lowered 
towards zero (see eg. This may seem surprising 

if one assumes that the mechanism allowing the mag- 
netisation to relax is thermal activation over a charac- 
teristic energy barrier AE. Namely, when T <C AE the 
relaxation rate should vanish since the Arrhenius factor 
for thermal activation, exp(— AE/T), goes to zero. The 
question then arises how relaxation can take place at a 
finite rate while the thermal activation factor exponen- 
tially approaches zero. 

A number of researchers have suggested that such a 
phenomenon is caused by quantum tunnelling of vortices 
through the barriers in the random pinning potential (for 
a theoretical review see [0). The above "quantum" ex- 
planation is very intriguing and in good agreement with 
some experimental results in compounds such as YBCO 
|,| or BSCCO f§. Also other even more exotic materi- 
als, such as organic superconductors j6), give good corre- 
spondences. However, in other non conventional systems 
such as heavy fermion superconductors, the theory of 
quantum creep is totally unable to describe the observed 
low T relaxation M. In fact, strong discrepancies are 
found in many other systems ranging from PCeCO crys- 
tals to YBCO/PBCO multilayers or YBCO and BSCCO 
films and crystals ||. One problem is || that the length 
of the tunnelling vortex segment L c (0) needed to fit the 
creep rate data, S, can be orders of magnitude larger that 
the one theoretically predicted by quantum creep theory 
j|. Also the experimental temperature dependence of 
the creep rate, S, is often very different from the one 
predicted by quantum theory (see ||). 

The above contradictory results suggests to look for 
additional descriptions of the anomalous low T magnetic 
relaxation. It is worth to stress that the observation of a 
non-vanishing constant creep rate in the limit T — > is 
found under very general circumstances: it does not cru- 
cially depend on the thickness of the sample ||] (i.e. on 
its dimensionality), nor on whether the pinning is caused 
by columnar defects or random point pins ||. Thus, the 
mechanism behind the low temperature creep seems to 



be of a fundamental and basic nature. 

We demonstrate below that also in a "classical" system 
(i.e., not "quantum") logarithmically slow glassy dynam- 
ics can naturally persist even at vanishing temperatures 
and can lead to the experimentally observed phenomenol- 
ogy. This is possible because the low T off-equilibrium 
dynamics consists of searching, among a very large num- 
ber, for a few "downhill" or "flat" directions in the free 
energy landscape. The number of these directions de- 
creases as relaxation proceeds though there always re- 
main some. They can be found only by collective co- 
operative rearrangements of the system, resulting in a 
slowing down of relaxation |10|. 

The model - We study a statistical mechanics model 
for vortex matter called a Restricted Occupancy Model 
(ROM) In the limit of zero temperature and infinite 
upper critical field, it reduces to a cellular automaton 
introduced in to study vortex avalanches. We use 
Monte Carlo dynamics which enabled the ROM to depict 
a unified picture of creep and transport phenomena in 
vortex physics, ranging from magnetisation loops with 
"anomalous" second peak, logarithmic relaxation, Bean 
profiles, to history dependent behaviours in vortex flow 
and I-V characteristics, to the reentrant nature of the 
equilibrium phase diagram JTH . The model also predicts 
the existence of a "glassy region" at low temperature with 
strong "aging" effects |[LlJ| . 

Here we use the ROM to study the magnetic relaxation 
rate, S, in the very low temperature limit. Interestingly, 
the model reproduces the experimental "anomalous" re- 
laxation and the observed behaviour of S @-||,|8| ■ 

The ROM model is described in full details in Ref. 
PH . A system of straight parallel vortex lines is coarse 
grained in the xy-plane by introducing a square grid of 
lattice spacing l of the order of the London penetration 
length, A |pH . The number of vortices on the i-th coarse 
grained plaquette is denoted by rij. The occupancy of 
each plaquette is a number larger than zero and, impor- 
tantly, smaller than N C 2 — B C 21q/4>q, where B C 2 is the 
upper critical magnetic field and <fto = hc/2e is the mag- 
netic flux quantum. The ROM model is thus defined by 
the following coarse grained vortex interaction Hamilto- 
nian Q: H = V ; / n,.\,,ii, - \ Y,i A a n i ~ J2i A i n i- 
The first two terms describe the repulsion between the 
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vortices and their self energy. On-site and nearest neigh- 
bour interactions are included, i.e., the interaction be- 
tween vortex lines with a separation greater than the 
London screening length (which by definition is close to 
lo) is ignored. We choose An = Aq = 1; A^ = A\ if 
i and j are nearest neighbours and Aij — otherwise. 
The last term in H describes a delta-distributed random 
pinning P(A*>) = (1 - p)S(AP) + pS{AP - A%). Interest- 
ingly, the general scenario of creep phenomena we de- 
scribe below does not depend on the details of pinning in 
the system (a fact in correspondence with experimental 
results In our model A sets the energy scale. Be- 
low we choose A = 1.0; A? = 0.3; N c2 = 27; p = 1/2; 
k* = Ai/A e [0.26,0.3]. 

The relaxation of the model is simulated by use of 
Monte Carlo dynamics on a square lattice in presence 
of a thermal bath of temperature T. The system is pe- 
riodic in the y-direction. The two edges parallel to the 
y-dircction arc in contact with a vortex reservoir. Par- 
ticles can enter and exit the system only through the 
reservoir, which plays the role of the external magnetic 
field. Hence the reservoir density, N ext , is used as the ex- 
ternal control parameter. We perform the following zero 
field cooled experiment: at a low temperature T we in- 
crease at constant rate 7 = ANq/t the reservoir density 
from zero to a working value N ext . We keep N ext fixed 
while we monitor the time dependence of the magneti- 
sation M — Ni n — N ex t. Here Ni n — Y) 4 n* / L 2 is the 
vortex density inside the system (of size L 2 Jl3[ ) . Time 
is measured in units of single attempted updates per de- 
grees of freedom of the lattice (see Ref. The data 
presented below are averaged over 128 realizations of the 
pinning background. 

In particular, we investigate the creep rate 
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as function of T, N ext and 7. In typical experiments 
the nature of the t dependence of M is such that S de- 
creases in time. So usually, one deals with an average 
creep rate, S a , in some given temporal window 
As shown in the upper inset of Fig.^J, in our model dy- 
namics M(t) at low temperatures behaves according to 
the known logarithmic interpolation formula (see Ref. 
0) found in experiments {§, namely: M(t) - M(0) ~ 
AM^l - [1 + f^ln(^)]- 1 ^ 1 }. Here, M(Q) is the 
magnetisation at the time of preparation of the sample 
(t = 0), AMjxj its overall variation, the exponent fi is 
consistent with 1, and to are fit parameters 0. 

Consistently, we define S a as the average value of S 
in the last time decade of our measures (i.e., for t G 
[10 5 , 10 6 ]). The present choice, analogous to those made 
for experimental data, is the most natural one and the 
general results presented below do not depend on it. 

For the reasons explained in the introduction, a very 
important physical quantity is the distribution, P{AE), 
of the energy barriers, AE, that vortices segments meet 



during their motion. Since at low T the system is typ- 
ically off-equilibrium, P(AE) is itself a (logarithmically 
slow) function of t and we consider its average over the 
last time decade of our measurements. 
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FIG. 1. Main frame The creep rate, S a , in the ROM 
model (k* = 0.28, 7 = 10 -3 ) for Next — 10 as a function of 
the temperature, T, in units of Aq. The error bars are of the 
size of symbols. The superimposed line is a linear fit. Inset 
The creep rate measured in BSCCO single crystal from Aupke 
et al. [5] for an external field of 880 Oe. 

Results - When the temperature is very low, the model 
exhibits the same kind of "anomalous" creep found in the 
experiments on superconductors. In Fig.|lj, we plot the 
creep rate, S a , as a function of T in a broad temperature 
range. For comparison we present equivalent experimen- 
tal measurements in a BSCCO single crystal (from Ref. 
II) as inset. The numerical values found for S a at low T 
in our model and in real samples are interestingly very 
similar. The temperature scales of the simulations and 
of real experiments can be compared by considering that 
the ratio T/Aq in our model is of the same order of magni- 
tude as T/T c in a real superconductor. This is seen from 
a comparison of the (T, N ext ) equilibrium phase diagram 
of our model with the equilibrium temperature-magnetic 
field, (T, H), phase diagram of, say, a BSCCO supercon- 
ductor (see Ref. @). 

In both experiments and simulations, S a approaches 
a finite value, S®, when T — > 0. In particular, we find 
that a linear fit of S a (T) in the low T regime is very 
satisfactory (see Fig.nl): 



S a (T)=S° a + aT 



(2) 



where both S® and a are a function of the applied field 
N e xt- We also note that in the present model S a (T) is non 
monotonous in T: in the higher T region it starts decreas- 
ing. This is also a known experimental fact p,^5|, we 
discuss it later on. The maximum in S a (T) is just above 
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a characteristic crossover "glassy" temperature, T g , de- 
fined in Ol. 
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very small decrease of S a with increasing 7 cannot be 
excluded: we show a fit to the form S° ~ Si + S2 111(7) > 
with 5i = 0.11 and s 2 = -2 • 10~ 4 . The fact that S° 
is pract ically independent on 7, far from being a proof 
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of the presence of equilibrium in the system, is due to 
the fact that at very low temperatures the characteris- 
tic equilibration time, T eq , is enormous (see Ref. 



So wheaever the driving rate, 7, is much larger than 
1 the off-equilibrium state and dynamics of the sys- 
teirh-aie essentially independent of 7. Stronger 7 effects 
have tp be expected when 7 gets closer to T~ q . In fact, 
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j^ftj^ (experimentally well known that at higher temper- 
atures the systems exhibit strong 7 dependent "mem- 
ory" effects [[ll 15 the signature of off-equilibrium 
dynamics. Actually, in the present model at low T, it is 
possible to show that r eq (T) diverges exponentially [jll| : 

• ferCO 7SJ ex P ("t 1 )- I n that region, the typical observa- 
^t?Mi time windows, t oos , are such that t O0S /T eq <C 1, and 
*tl™jjjSy.§tem is in the early stage of its off-equilibrium re- 
laxation from its initial state. This is schematically the 
origin ol f the flattening of S a at very low T. Notice that, 
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FIG. 2. Main frame The zero temperature limit of the 
creep rate, S a , in the ROM model as a function of the applied 
field N ext (for T = 1(T 4 and 7 = 10~ 3 ). The superimposed 
curve is a power law to guide the eye. Inset top: The relax- 
ation of the magnetisation, M(t), in the model («* = 0.28, 
7 = 10 -3 ) for N ex t = 10 and T = 0.25 as a function of time. 
The continuous line is the logarithmic fit of the text. Inset 
bottom: S a is here plotted as a function of the sweep rate, 



7, of the external field for N ex t 
weak dependence is found. 



10 and T = 10 . A very 



In our model by varying the applied field we find a 
range of values for S® very similar to experimental ones 
[Q-Q (see Fig.||). In particular, 5° seems to decrease 
on average by increasing the field N ext . The overall be- 
haviour can be approximately interpolated with a power 
law: S° a (N ext ) ~ (N ext /N )-\ where, for k* = 0.28, 
Nq — 0.01 and x ~ 0.6. As shown in Fig.||, the presence 
of a small exponent x implies that sensible variations in 
Sq can be seen only by changing N ext of orders of magni- 
tude. Note that the dips in the S a (0) versus N ext data in 
Fig. H at certain values of N ext (namely around 3, 13, and 
18) are statistically significant. They are related to the 
low field order-disorder transition, the 2nd peak transi- 
tion and the reentrant high field order-disorder transition 
respectively ]TT| . 

In the lower inset of Fig.|^, we show that S a is essen- 
tially independent of the ramping rate, 7 (the values 
shown are for N ext = 10 and T = 10~ 4 ). This is an 
other typical experimental observation M. However, a 



if one could observe the system for an exponentially long 
?ime, i.e., if t (, s /T eq » 1, then the creep rate, S a , would 
indeed go to zero. 

The above scenario is clarified by the analysis of the en- 
ergy barrier distribution function, P(AE), recorded dur- 
ing the system evolution at very low T. Such a quantity 
also clearly shows the simple mechanical origin of the 
anomalous creep found at very low temperature in the 
present model. The function P(AE) (where AE is in 
units of Aq), recorded at T — 10~ 4 , is plotted in Fig.|| 
for two values of the applied field, N ext . We always find 
that P(AE) has support also on the negative axis. This 
is the mark of the off-equilibrium nature of the evolution 
on the observed time scales. The presence of a P(AE) 
which extends down to negative values also explains the 
presence of the recorded relaxation at low T: in the con- 
figuration space the system can still find directions where 
no positive barriers have to be crossed. The insert in Fig. 
|3| clarifies the mechanism behind the relaxation. Here, we 
plot the signal A(t) defined, for each single Monte Carlo 
(MC) step t, in the following way: A(t) = if the MC 
trial is rejected; A(t) — 1 if the trial is accepted and 
the energy reduced, i.e. AE < 0; and finally A(t) = 2 
when a trial is accepted with AE > 0. We plot two se- 
quences of trials. One for < t < 500 was measured 
at the early stage of the relaxation, the second sequence, 
placed at the interval 500 < t < 1000, is measured at 
the late stage of the relaxation. Most trials are rejected 
(A(t) = 0) and only once in a while the system does 
manage to find a route pointing downhill in the energy 
landscape. A(t) = 2 never occurs. As time proceeds 
fewer and fewer "negative channels" are available to the 
relaxation and a decrease in the density of the spikes in 
A(t) is observed. 

As the temperature is increased thermal activation 
over positive energy barriers will become possible as the 
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Arrhenius factor exp(— AE/T) assumes a non-vanishing 
value for an appreciable range of barrier values AE < T. 
When this happens the relaxation will occur sufficiently 
fast to allow one, within the experimental time window, 
to closely approach the equilibrium configurations where 
the vortex density profile is more or less flat and relax- 
ation ceases, hence S a goes down, as seen in experiments 
L5| and in Fig. 1. 
-1 
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FIG. 3. The energy barrier distribution, P(AE), recorded 
at T = 10~ 4 for N ext = 25,4 (resp. circles,diamonds) . No- 
tice that P(AE) extends down to negative values. Inset We 
plot the function A(t) defined in the text to monitor the ac- 
tivity during relaxation. In the interval [0, 500] we show A(t) 
recorded at the beginning of the run and in [500, 1000] A(t) 
in the last steps of the same run. 

Finally, we stress that slow off-equilibrium relaxations 
at very low temperatures are also observed in glass form- 
ing liquids |10[j . In that cases too, no activation over bar- 
riers occurs and the system simply wanders in its very 
high dimensional phase space through the few channels 
where no energy increase is required. 

Discussion - We have above demonstrated that the 
phenomenological behaviour of the creep rate at low 
temperatures can be understood in terms of the off- 
equilibrium nature of the inhomogeneous vortex density 
profile produced in magnetic creep experiments. In fact, 
cooperative mechanical rearrangements, possible even at 
very low T (where thermal activation over positive bar- 
riers can be negligible), dominate the phenomenon [fl0|| . 
In this perspective, it is very important to stress that 
the system's equilibration time at very low temperature 
is much larger than any experimentally accessible time 
window p| . Accordingly, we can say that experimen- 
tal findings do not enforce the interpretation in terms of 
macroscopic quantum tunnelling of vortices. 

Relaxation due to quantum tunnelling might be 
present along with the mechanical relaxation discussed 



above. It is then important to ask how compelling the 
quantum tunnelling interpretation is. The theory of 
quantum tunnelling assumes a London picture and treats 
the position of the vortex core as the variable that is 
able to tunnel. It is not entirely clear if this is the right 
level at which to introduce quantum fluctuations, but, 
more importantly, the quantum tunnelling description 
assume the existence of a characteristic energy barrier 
, not a time dependent distribution of barrier heights 
as typically found in many-body systems relaxing off- 
equilibrium. Finally, the quantum tunnelling description 
also tacitly assumes the existence of a static equilibrium 
state in which barriers are always positive. As we have 
clearly shown above, this is typically not the case and a 
dynamical approach is more appropriate. 

The present scenario, where off-equilibrium phenom- 
ena dominate the anomalous low T creep, could be exper- 
imentally verified by the discovery of "aging" phenomena 
|i~Tf like those recently observed in Jl7],[l8| . 
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